AD-AO01  667 

UNCLASSIFIED 


SCIENCE  APPLICATIONS  INC  MCLEAN  VA  F/0  20/1 

DEPTH  DEPENDENCE  AND  VERTICAL  DIRECTIONALITY  OF  AMBIENT  SEA  NOI— ETC 
AUG  79  M  W  RENNER »  R  C  CAVANAGH  N00014-79-C-0145 

SAI-81-151-WA  NL 


j>EPTH  DEPENDENCE  AND  JERTICAL 
f  DIRECTIONALITY  OF  ^MBIENT  JSEA  JfOISE 


SAI-81-151-WA 


Prepared  by: 

William  W. h  enner 
Raymond  C.J  Cavanagh 


_ _ _ _ _  Prepared  lor: 

•■'f '  •  (.  s't-i'yi » *J 

Oeean  ^og^jms  Of f ice  (Code  500) 


Naval  Ocean  Research 
Development  Activity 
NSTL,  Mississippi 


SCIENCE  APPLICATIONS,  INC. 


8400  Westpark  Dr. 

McLean,  Virginia  22101 
(703)  821-4300 


v.pnw<l  lot  yubllc  \ 

"Distribution 


TABLE  OF  CONTENTS 


Page 

Section 

1 

INTRODUCTION  AND  BACKGROUND 

1-1 

1.1 

Introduction 

1-1 

1.2 

Background 

1-2 

1.3 

Outline 

1-4 

Section 

2 

NOISE  MODEL  LIMITATIONS  AND 
APPROXIMATIONS 

2-1 

2.1 

The  Model  and  Assumptions 

2-1 

2.2 

Omnidirectional  Noise  Versus  Depth 

2-5 

2.3 

Surface-Image  Interference  and 
Diffraction  Effects 

2-8 

2.4 

Application  to  Transmission  Loss 

2-11 

Section 

3 

EXAMPLES  OF  DEPTH  DEPENDENCE  FOR 
CANONICAL  DIRECTIONALITY 

3-1 

3.1 

Isotropic  Vertical  Directionality 

3-1 

3.2 

Vertical  Directionality  as  a 

Step  Function 

3-2 

3.3 

Vertical  Directionality  with  a 
Singularity 

3-3 

Section 

4 

INVERSE  PROBLEM 

4-1 

4.1 

Analytic  Solution 

4-1 

4.2 

Error  Minimization  Algorithms 

4-5 

4.2.1 

Least-Squares  Estimate 

4-9 

4.2.2 

Linear  Program 

4-10 

Section 

5 

EXAMPLES  OF  EFFECTS  OF  NOTCH 

FILLING  ON  DEPTH  DEPENDENCE 

5-1 

Section 

6 

SUMMARY  AND  CONCLUSIONS 

6-1 

1 

1 

i 

i 


i 


Section  1 


INTRODUCTION  AND  BACKGROUND 


1.1  Introduction 


The  performance  of  passive,  low-frequency  acoustic 
sensors  is  largely  limited  by  prevailing  ambient  sea  noise. 
In  the  past,  the  Navy  has  focused  research  on  the  omni¬ 
directional  (omni)  spectrum  levels  and  depth  dependence  of 
the  noise  field  to  optimize  design  and  deployment  of  shallow 
and  mid-water  systems  (e.g.,  sonobuoys).  Because  of  recent 
developments  in  sonar  arrays  with  vertical  apertures, 
there  is  now  considerable  interest  in  the  vertical  direc¬ 
tionality  of  ship  and  wind-generated  noise.  This  report 
develops  and  applies  a  theoretical  relationship  connecting 
these  two  important  properties  of  sea  noise:  depth  depend¬ 
ence  of  the  omni  level  and  vertical  directionality  at  a 
single  depth. 

The  relationship  is  based  on  an  analytical  result 
which  was  documented  a  few  years  ago  (Reference  1-1),  but 
which  has  apparently  not  been  exploited  to  date.  It  is 
somewhat  limited  by  assumptions: 

the  noise  field  is  an  "average"  in  time, 

-  geometric  acoustics  applies, 

the  depth  regime  is.  limited  to  a  single  "sound 
channel"  (see  Subsection  2.1). 
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However,  the  result  Is  not  limited  to  range-independent 
environments  and  has  applications  to  many ' problems , 
including: 


analysis  and  interpretation  of  noise  data, 

insight  into  anomalies  in  the  noise  field 
(e. g. ,  the  "notch") , 

optimization  of  computer  noise-prediction 
models. 

The  approach  applies  as  well  to  range-averaged 
transmission  loss  (see  Subsection  2-4).  Specifically,  under 
the  assumptions  listed,  the  relationship  between  range- 
averaged  transmission  loss  as  a  function  of  depth  and 
vertical  arrival  structure  is  similar  to  the  relationship 
between  the  corresponding  properties  of  sea  noise  discussed 
in  this  paper. 

The  present  work  was  motivated  by  an  effort 
(Reference  1-2)  to  improve  the  noise  module  of  a  verti¬ 
cal-array  performance-prediction  model  used  for  Fleet 
support.  The  result  allows  omnidirectional  depth-dependence 
data  to  give  clues  about  the  vertical  directionality  of  the 
noise  field.  More  detail  on  this  background  follows. 

1.2  Background 

A  new  generation  of  passive  sonobuoys  which  employ 
vertical  line  arrays  (VLA)  can  achieve  significant  gains  in 
signal-to-noise  ratio  (compared  to  an  omnidirectional 
sensor)  by  discriminating  in  favor  of  specific  signal 
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arrival  angles  and  against  ambient  noise  arriving  from  other 
directions.  A  performance  model  for  such  a  sensor  requires 
a  prediction  of  the  vertical  directionality  of  the  ambient 
noise  at  the  sonobuoy  depth.  The  current  Fleet-support 
model  (ASRAP)  estimates  vertical  directionality  with  an 
approach  due  to  Talham  (Reference  1-3)  and  extended  from 
FANM  (Reference  1-1).  Details  about  the  noise  component 
of  ASRAP  can  be  found  in  Reference  1-4. 

A  problem  with  the  current  ASRAP  noise  module  is 
the  chronic  prediction  of  a  "notch"  (absence  of  energy)  at 
angles  near  the  horizontal.  Because  of  the  assumption  of  a 
range-independent  environment  (i.e.,  sound  speed  and  bathy¬ 
metry  do  not  change  with  range)  and  the  use  of  geometric 
acoustics  (i.e.,  no  diffraction  or  scatter),  the  model  will 
predict  a  null  in  the  noise  directionality  whenever  the 
sound  speed  at  the  receiver  depth  is  less  than  the  sound 
speed  at  some  point  above  it.  In  contrast,  models  which 
account  for  range-dependence  indicate  that  the  notch  is 
often  partially  filled  and  in  many  cases  significantly 
filled.*  Furthermore  this  is  qualitatively  supported 
by  several  noise  vertical  directionality  measurements 
(Refs.  1-5  and  1-6). 

Because  of  the  limited  quantity  and  quality  of 
noise  directionality  data,  additional  information  was  sought 


*  Diffraction  and  scatter  are  estimated  to  be  of  secondary 
importance  at  frequencies  above  25  Hz. 
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to  aid  in  determining  the  extent  to  which  the  notch  is 
filled.  The  relationship  discussed  in  this  report  allows 
ambient-noise  depth-dependence  measurements  to  be  applied 
to  the  problem,  and  this  report  is  chiefly  concerned  with 
how  such  data  can  be  used  to  estimate  vertical  direction¬ 
ality. 


1.3  Outline 


The  relationship  between  noise  depth  dependence 
and  vertical  directionality  is  considered  in  both  direc¬ 
tions.  Section  2  discusses  the  assumptions  and  limitations 
of  the  approach,  and  then  develops  formulas  and  approxi¬ 
mations  for  determining  the  depth-dependence  from  the 
vertical  directionality  at  a  specific  depth.  Section  3  then 
derives  the  depth-dependence  function  for  several  canonical 
directionality  functions,  and  Section  5  provides  examples  of 
one  of  these  depth-dependence  functions  for  a  variety  of 
environments. 

Section  4  deals  with  the  inverse  problem:  de¬ 
termining  the  vertical  directionality  from  the  depth  depend¬ 
ence.  Two  different  approaches  are  considered:  an  analytic 
solution  is  developed  for  the  case  that  the  depth  dependence 
has  a  particular  representation;  an  error-minimization 
algorithm  is  applied  when  the  vertical  directionality  has  a 
particular  form. 
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Section  2 


NOISE  MODEL,  LIMITATIONS  AND  APPROXIMATIONS 

This  section  gives  relationships  between  omni¬ 
directional  ambient  noise  and  vertical  directionality,  as 
functions  of  depth,  based  on  a  rather  general  model.  The 
limitations  and  assumptions  associated  with  the  model 
are  described,  and  approximate  formulas  derived. 

2 . 1  The  Model  and  Assumptions 

Begin  by  assuming  that  the  directionality  of 
the  noise  field  Ns(0,<|>;zo)  is  known.  Here,  N  represents 
equivalent-plane-wave  intensity,  per  steradian,  at  depth  zQ, 
for  vertical  angle  0  (measured  from  the  horizontal)  at 
azimuthal  angle  <J),  in  a  small  band  of  frequencies.  Next 
suppose  that  geometric  acoustics  gives  an  accurate  approxi¬ 
mation  of  the  field  and  that  the  medium  is  range-independent 
in  a  neighborhood  of  the  receiver  location  (specifically,  in 
a  region  including  a  ray's  range-cycle;  say,  35  miles  at 
most).  Azimuthal  homogeneity  is  not  required. 

Under  such  assumptions,  the  vertical  direction¬ 
ality  (per  steradian)  can  be  viewed  as  the  average  in 
azimuth  of  arrivals; 

2tt 

Ns<6-V  -  J  Ms«M:v  d*.  (2-D 
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Furthermore,  since  Snell's  law  holds  locally,  the  noise 
reaching  the  receiver-depth  z0  at  angle  9  will  reach 

A 

depth  z  at  angle  6 , 


8 


sgn(0 ) *  arccos 


c(z) 

c<zo> 


cos 


(2-2) 


where  c(z)  is  the  sound  speed  at  depth  z. 
then  the  path  will  not  reach  depth  z] . 


llf  Ifv  cos  el>1’ 


Continue  by  assuming  that  the  sources  of  noise 
are  distributed  as  a  continuum,  homogeneous  in  range 
(radially  from  the  receiver)  over  a  ray's  range  cycle. 
For  wind  sources,  this  is  reasonable.  For  ship  sources, 
view  the  noise  field  as  an  average  in  time  (and  hence  over 
source  or  receiver  location).  This  additional  premise 
allows  arrivals  to  be  translated  in  range,  with  the  only 
error  resulting  from  volume  attenuation  for  the  difference 
in  path  lengths  (negligible  at  low  frequency).  Such  trans¬ 
lations,  illustrated  in  Figure  2-1,  lead  to  the  important 
relationship, 


ns(9,$;z)  <*  Ns(e,$;z0) 

and  for  vertical  directionality  (applying  (2-1)); 

Ns(e;zo).=  Ns(6A;z) 

where  6  and  §  are  related  by  (2-2). 


(2-3) 


(2-4) 
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The  consequences  of  these  observations  are  sig¬ 
nificant.  For,  once  ^s(B  ,  <p  ;z0)  is  known  for  all  0  at  a 
given  depth  z0,  it  can  be  estimated  at  any  depth  z,  provided 
that  every  ray  path  that  reaches  z  also  reaches  zQ.  For 
the  case  that  z0  is  a  local  minimum  in  the  sound  speed 
profile,  this  requirement  is  met  if  z  is  in  the  same  "sound 
channel"  as  zQ,  i.e.  there  is  no  local  maximum  between 
z  and  zQ.  Furthermore  if  the  environment  is  assumed  to  be 
range  independent,  this  restriction  can  be  lessened  to:  no 
local  maximum  between  z  and  z0  that  is  greater  than  the 
sound  speed  at  the  source  depth  (usually  the  surface).  It 
is  important  to  note  that  while  all  rays  that  reach  depth  z 
also  reach  zQ  under  these  conditions,  the  converse  is  not 
true.  This  can  be  seen  in  that  the  transformation  from 
is  defined  only  on  the  domain  a  (z) < j 0 |<J  where 


As  will  be  observed  time  and  again  in  what 
follows,  under  the  assumptions  of  the  model  the  depth 
dependence  of  noise  directionality  is  directly  related  to 
depth  changes  in  sound  speed.  In  fact,  if  zj  and  Z2 
are  in  the  same  "sound  channel"  and  if  c(zj)=c(z2),  then 
the  noise  directionality  (N(0;z))  is  the  same  at  zj  and 
Z2-  Thu§,,  to  get>  new  information  about  the  noise  field, 
samples  must  be  taken  at  depths  with  different  sound 
speeds.  To  g-o  one  step  further,  a  case  of  constant  sound 


*  Note  that  o<c(z0)/c(z)  <  1_.  a(z)  is  taken  as  the  value 
of  the  arccos  between  0  and  tt/2. 


speed  should  be  viewed  as  degenerate;  for  then  the  direc¬ 
tionality  and  omni  levels  are  also  constant  with  depth  -  so 
that  directionality  cannot  be  determined  from  depth 
dependence . 

In  summary,  the  noise  vertical  directionality 
depends  on  depth  in  a  simple  way  (equations  2-3  and  2-4) 
under  the  assumptions  that: 

Plane-wave  geometric  acoustics  applies. 

Surface-image  interference  and  diffraction 
effects  are  ignored  (but  see  Subsection  2.3). 

Volume  attenuation  is  small  over  15  miles 
(low  frequencies). 

-  The  environment  is  range-independent  within 
about  30  miles  of  the  receiver  location. 

Sources  of  noise  are  continuously  distributed, 
and  approximately  homogeneous  over  ranges  of 
30  miles. 

-  The  formulas  are  applied  for  depths  within 
the  same  "sound  channel". 

Note  that  the  -assumptions  about  the  environment  and  sources 
are  equivalent  to  those  of  FANM  II  (Reference  2-1),  and  are 
consistent  with  any  noise  model  which  uses  range-averaged 
transmission  loss  in  a  range-dependent  environment  (e.g., 
ASTRAL  in  an  adiabatic  case). 
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2.2  Omnidirectional  Noise  Versus  Depth 


Under  the  additional  assumption  that  arrivals 
add  incoherently,  the  omnidirectional  noise  at  depth  z0 
can  be  calculated  directly  as: 

N(z0)  =  JJh s  (0,4>;zo)  d (2 

UNIT' 

SPHERE 


N  C  ©  ,  4> ;  z  )  cos  0  d0  d<}>. 
s  o 


(2-5) 


In  terms  of  vertical  directionality  (see  (2-1)), 

r/2 

Ns(0;zo)  ^os  0  d0 . 

-tt/2 


N(zq)  *  2tt 


•/-7T  J 


(2-6) 


(Again,  notice  that  if  c(z)  is  constant,  then  so  are  Ns(0;z) 
and  N(z)  as  functions  of  z). 


When  c(zQ)  is  a  local  minimum  in  the  sound  speed 
profile  and  z  is  in  the  same  sound  channel,  then  the  omni¬ 
directional  depth  dependence  can  be  calculated  from: 


N(z) 


N  (6;  z)  cos  0  d0 


^  a 

z«)  cos  0  <*0 
o  O 


(2-7) 


2-6 


where 


FCe.z) 


sin6  cos6 


cos  [ot C z )]  ^cos2  [  a ( z )]  -  cos2e 


(2-8) 


and 


a(z) 


arccos 


C(2p> 

C(Z) 


In  the  case  that  Ns(e;z0)  is  symmetric  about  the  horizontal, 
(2-7)  simplifies  to 

N(2)  '  “”/a(z)  <2-9 

Examples  of  the  application  and  interpretation 
of  these  results  are  found  in  Section  3. 
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Surface-Image  Interference  and  Diffraction 
Effects 


To  this  point,  the  changes  in  vertical  direction¬ 
ality  with  receiver  depth  have  been  attributed  completely  to 
the  geometric  change  in  ray-path  structure,  calculated 
directly  from  Snell's  law  in  a  neighborhood  of  the  re¬ 
ceiver.  Moreover,  in  the  last  subsection  a  formula  for  the 
depth  dependence  of  the  total  (omnidirectional)  level  was 
derived  under  the  additional  assumption  that  the  noise 
consists  of  plane-wave  arrivals  which  add  on  a  random-phase 
basis  (incoherently).  The  following  discussion  is  concerned 
with  two  mechanisms  which  can  modify  the  results  obtained 
from  the  simplified  ray  model:  surface-image  interference 
(SII)  and  diffraction. 

Concentrate  first  on  noise  arrivals  from  a  single, 
distant  point  source.  If  these  arrivals  were  added  in 
phase,  a  complex  multipath-interference  field  would  be 
observed  as  the  source  or  receiver  moved  in  range  or  depth. 
Because  of  the  assumptions  of  time-  or  range-averaged  noise 
(as  reflected  in  the  noise-source  distribution)  associated 
with  the  model  of  this  paper,  only  those  fluctuations  which 
persist  in  range  need  be  considered;  the  others  are  smoothed 
in  the  averaging  process.  Following  the  rationale  of 
FACT  (Ref  2-2),  we‘ limit  the  candidate  interfering  paths  to 
pairs  of  long-range  arrivals  which  differ  in  their  history 
by  a  single  surface  reflection  near  the  receiver.  The 
resulting  effect  can  substantially  increase  low-frequency 
transmission  loss  for  receivers  near  the  surface,  and  is 
usually  called  "Surface-Image  Interference"  (in  special 
cases,  "Lloyd's  Mirror"). 
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To  account  for  the  SII  Mechanism,  first  note 
that  the  noise  directionality  is  unaltered;  phases  can 
be  assigned  to  the  arrivals,  but  until  the  arrivals  are 
summed  there  is  no  manifestation  of  the  interference 
(There  is  an  implicit  assumption  here  that  arrivals  from  two 
"sources"  add  incoherently).  Focus  than  on  the  calculation 
of  the  omni  level  (N(z)).  At  least  two  approaches  make 
sense: 

•  Arrivals  from  above  and  below  with  the 
same  (absolute-value)  angle  are  presumed 
to  interfere.  The  well-known  approximation 
can  be  used  to  modify  the  sum: 

I2  -  4Il  sin2  ,  (2-10) 

where  Ij  is  the  intensity  for  a  single  arrival 
at  angle  8 ,  w  is  radial  frequency,  c  is  sound 
speed,  and  z  is  receiver  depth.  More  precise 
formulas  including  the  influence  of  the  local 
sound-speed  variations  with  depth  are  avail¬ 
able  (see,  e.g.,  Ref.  2-3).  Vhen  the  changes 
of  I2  with  depth  become  rapid  (as  w ,  z,  or 
0  increase),  I2  is  replaced  by  the  incoherent 
(RMS)  sum. 

•  The  TAPPS  noise  module  (Ref.  2-4)  predicts 
noise  depth  dependence  for  towed  array 
systems  and  attributes  all  variation  to  SII. 
Since  the  array  is  presumed  shallow,  the 
frequency  low,  and  arrivals  from  ships  at 
small  angles,  such  a  treatment  is  not  unrea¬ 
sonable.  Note,  however,  that  the  sound  speed 


profile  at  the  array  is  used  only  to  calculate 
limiting  ray  arrivals  from  the  surface  and 
bottom. 

TAPPS  considers  two  regimes  of  noise  arrivals: 
paths  which  do  not  suffer  bottom  attenuation, 
and  those  which  do.  In  each  case,  an  average 
of  the  sum  of  (2-10)  over  the  appropriate 
angular  aperture  is  used  to  find-  omni  noise 
depth  dependence.  Only  in  tbe  bottom-bounce 
case  are  the  individual  arrivals  within 
the  aperture  presumed  to  change  in  intensity 
with  angle.  The  TAPPS  algorithm  is  nearly 
the  same  as  that  of  tbe  FACT  shallow-water 
routine  (Ref.  2-1). 

Either  of  these  approaches  can  be  applied  directly  to 
the  model  described  in  Subsection  2.2  to  account  for  SII. 
It  should  be  emphasized  that  the  effect  is  important  only 

for  small  values  of  aiz  sin  8  . 

c 

Other  physical  mechanisms  not  included  in  the. 
model  are  those  related  to  "wave"  or  diffraction  phenomena. 
They  can  be  important  at  low  frequencies,  but  are  not 
predicted  by  geometric  acoustics.  For  example: 

•  The  noise  field  actually  extends  into  "shadow" 
regions  bounding  the  "sound  channel."  Tbe 
intensity  and  extent  of  the  diffracted  field 
are  inversely  proportional  to  frequency 
and  sound-speed  gradient. 


*  • 


•  Ray  arrivals  which  cycle  within  a  constrained 
depth  interval  (e.g.,  in  a  surface  duct  or 
about  a  channel  axis)  are  subject  to  model 
attenuation  as  frequency  decreases  (low-fre¬ 
quency  cut  off). 

Such  effects  can  be  included  in  a  geometric-acous¬ 
tic  framework,  as  has  been  done  in  the  ASTRAL  (Ref.  2-5) 
and  other  transmission-loss  models. 

In  summary,  SII  and  diffraction  can  significantly 
modify  noise  depth  dependence  at  low  frequency.  There  are 
straightforward  methods  for  estimating  the  contribution  of 
each.  The  remainder  of  this  paper  will  concentrate  on  the 
geometric  aspects  introduced  in  Subsection  2.1  without 
further  consideration  of  these  other  phenomena. 


2.4  Application  to  Transmission  Loss 


I 


The  results  of  this  paper  apply  almost  without 
change  to  range-averaged  transmission  loss.  Specifically, 
suppose  S  (6;z^)  represents  transmission  arrival  structure 
at  the  receiver  depth  z^,  averaged  over  an  interval  of 
source  range  exceeding  the  longest  ray  period  (say,  35 
nm) .  Make  assumptions  similar  to  those  listed  in  Sub¬ 
section  2.1: 


•  Plane-wave,  geometric  acoustics  applies. 

•  SII  and  diffraction  are  ignored. 
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•  Volume  attenuation  is  small  over  15  nm. 


•  The  "source"  or  receiver  is  distributed 
over  35  nm  (i.e.,  TL  is  averaged  over  that 
range) . 

•  Formulas,  are  applied  for  depths  at  the 
receiver  within  the  same  "sound  channel". 

•  The  environment  is  range-independent  within 
about  35  nm  of  the  receiver. 

In  that  case,  formula  (2.4)  applies  to  : 


for 


S(6;zo)  =  S (8;z) 


8  =  sgn( 8 ) 


arccos 


(c(z^)  cos6)' 


Furthermore,  the  "total"  (as  received  by  an  omnidirectional 
receiver)  sound  intensity  S  (z)  satisfies  the  analogy  of 
(2-6): 


A 

zo)  d0. 


The  only  real  difference  between  the  transmission  and 
noise  problem  is  in  the  dimension  of  the  sound  field. 


The  formulas  can  be  applied  as  in  the  case  of 
noise  to  relate  arrival  structure  and  depth  dependence. 
The  last  of  the  listed  assumptions  should  be  emphasized: 
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it  precludes  application  of  the  result  to  receivers  located 
in  locally  range-dependent  environments  (e.g.,  it  could  not 
be  applied  to  the  ASTRAL  problem  with  receiver  over  a 
sloping  bottom). 

Finally,  there  is  also  an  obvious  application 
to  TL  and  AN  computer  models,  already  known  but  not  always 
recognized.  In  practice,  rays  are  usually  traced  from 
receiver  to  source.  If  the  receiver  is  placed  at  the 
channel  axis,  and  the  individual  ray  intensities  stored, 
then  the  formulas  for  angle  transformation  allow  the  arrival 
structure  or  directionality  to  be  constructed  at  other 
depths  within  the  same  "sound  channel"  without  additional 
ray  tracing.  The  only  error  suffered  is  in  volume  atten¬ 
uation  over  at  most  one  ray  cycle. 


Section  3 


EXAMPLES  OF  DEPTH  DEPENDENCE  FOR 
CANONICAL  DIRECTIONALITY 


For  certain  realistic,  albeit  simplistic,  vertical 
directionality  functions,  the  noise  depth-dependence  func¬ 
tion  can  be  found  analytically.  This  section  .derives  the 
depth  dependence  functions  for  three  such  cases.  In  all 
cases  c(z0)  is  a  local  minimum  in  the  sound  speed  profile 
and  z  is  restricted  to  the  same  sound  channel  as  zQ  (i.e., 
no  local  maximum  between  zQ  and  z). 

3.1  Isotropic  Vertical  Directionality 

For  frequencies  at  which  the  wind  noise  dominates, 
some  measurements  show  little  vertical  directionality  in 
the  ambient  noise.  Hence  consider  the  case  where, 


Ns(6;Zo)  =  N,  -tt/2  <  0<tt/2. 


Then  the  omnidirectional  noise  intensity  at  zQ  is. 


N(z  )  =  2 it  /^/2  N  (■© ;  z  )  cos6  d0 
°  •'-tt/2  s  ° 

=  4ttN, 
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and  at  z, 


N(z)  =  4  it 


f  tt/2 

J  ot(z) 


Ns(6;  zo).F(0,z)d8 


*7T  /  2 

=  4-rrN#  F(0,z)d0 


£/2 
.Z) 


=  4ttN 

^os^ 

a(z) 

2 

-COS  0 

cos 

a(z)l  . 

tt/2 

0=a(z) 


4ttN. 


Thus  an  isotropic  vertical  noise  pattern  yields  a  constant 
depth-dependence  function. 

3.2  Vertical  Directionality  as  a  Step  Function 

This  example  illustrates  the  "notch"  problem 
mentioned  earlier.  Let  cs  be  the  sound  speed  at  the  sur¬ 
face  and  cB  be  the  sound  speed  at  the  bottom  with 


cC2o)<  cs  -  CB 


0g  =  arccos 


c(zQ) 


0g  =  arccos 


c(zQ) 


L  cb  J 


(3-2) 
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Noise  models  which  place  the  source  at  the  surface  and 
presume  range-independent  environments  will  predict  that 
Ns(0 ,Zo)  is  greatest  when  6S£|0|<  0g  and  near  zero  when 
0<_  I9!*11  6S.  Consider  then  the  case  (see  Fig.  3-1)  for 
which, 

N1,o<|eJ<es 

Ns(0.zo>  =  N2*es-ie'l<6B 

N3>0B<|e|<TT/2  . 


Then, 


N(zq)  =  2  it 


N  (0;z  )  cos  0d0 
s  o 
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IT/ 

— 


-e, 


-e„ 


0 


0, 


B  'S  S  'B 

Vertical  Angle  at  the  Receiver 


n/ 

'o 


Figure  3-1 

Vertical  Directionality  as  a  Step  Function 
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.  J 


Vertical  Directionality  with  a  Singularity 


e, 


=  arccos 


Let  c5  denote  the  sound  speed  at  the  surface, 
/ c(zn)  \  r  Co 

:OS  \~ — /  ’  and  0G(6)=  arccos 

>e  models  (e.g.,  FANM  or  ASRAP  Refen 


cos  0 


In 


Talham-type  models  (e.g.,  FANM  or  ASRAP  References  1-1  and 
1-4),  an  isotropic  source  directivity  results  in: 


*  -ii5g:(9T  as  6-6s 


so  that 


lim  N  (0;z  )  - 
6-6  g 


Now 


COS0 


sine. 


COS  ©  -  COS2©  , 

s 


so  consider 


the  noise  directionality  function  (Fig.  3-2) 
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Ambient 


•  Angles  at  the  Receiver 

Figure  3-2 

Vertical  Directionality  with  a  Singularity 
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Then 


“  N1sin8s  . 


reB 

N2coses  I _ sin 

coso(z)  |  <003*8^008*0) 


sin8cos9d8 _ 

*8)^(cos*a(z)-cos*8)^ 


+  N3(l-sin6s) 


where  0  and  0R  are  as  defined  in  (3-2).. 


Letting  x  **  cos  ^9.  the  integral  above  becoaes 


'COS20. 


—  - -  ^ 

(COS20  -x)i(cos2a(z)-x)i 


'cos20. 


[1  I OOS*  0B 

2(cos*es-x)’(cos*o(z)-x)i  +  2x  -  cos*0s  -  cos2a(z)  I 
,  -*lx=*cos*e 


[•  cos*8s-cos*aCz; 

2(cos*0s-cos28B)“(cos,a(z)-cosJ 0B)^  + 


2cos20o-cos*0  -cos*o(z) 
p  s 


(3-4) 


*See  Reference  3-1 


Now  for  c(z)  f  cs,  N(z)  is  finite.  However  as  c(z)  — 
cs  then  cosa(z)  —  cos0s  and  thus  the  expressibn  in  (3-4) 
approaches  infinity.  The  reason  for  the  singularity  is 
discussed  in  Reference  1-1;  it  is  a  result  of  the  fact 
that  range-averaged  transmission  loss  calculated  with  ray 
acoustics,  has  a  finite  average  at  smooth  caustics  but  not 
at  horizontal,  cusped  caustics.  The  latter  occur  when 
source  and  receiver  are  at  depths  with  the  same  sound 
speed . 
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Section  4 


INVERSE  PROBLEM 

To  this  point  the  discussion  has  focused  on 
how  the  ambient  noise  depth  dependence  can  be  approximated 
from  the  noise  directionality  at  a  fixed  depth.  Of  equal 
importance  is  the  inverse  problem:  When  can  the  direction¬ 
ality  be  determined  from  the  depth  dependence?  This  section 
deals  with  two  approaches,  the  first  an  analytic  solution 
for  depth-dependence  functions  of  a  particular  form  and  the 
second  an  error  minimization  algorithm  for  directionality 
functions  similar  to  those  of  subsection  3.2. 


Analytic  Solution 


CASE  1  Symmetry  in  8 


Assume  Ns(0,zo)  is  symmetric  in  e  for  zQ  fixed 
at  the  axis  (minimum  sound  speed)  of  a  "sound  channel." 
Then  recalling  (2-9) 


N(2)  =  4n  I  Ns( 9  »  zq ) 


ra(z) 


cosa(z) 


sin0  cosB 


/cos2a(z)-cos2e 


(4-1 


where  a(z)  =  arccos 


/cM  . 

\  C(z)/'wl 


with  0  <_  ct(z)  <  tt/2, 


This  relationship  can  be  viewed  as  a  mapping 

of  vertical  directionality  functions  {Ng(8;  zQ)}into  omni 

depth  functions  (N(z)).  An  analytic  solution  to  the  inverse 

mapping  (from  (N(z)}  to  { N  ( 9 ;  z_)})  is  sought. 

s  o 
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Recall  (Section  2)  that  whenever  c(zj)  -  c(z2> 
and  zj,Z2  are  in  the  same  "sound  channel",  then 

Ng(e ;  z1)  =  Ns(0;  z£) ,  for  all  9, 


and 

N(z1)  =  N(z2). 


It  makes  sense  then  to  define  c“*  as  a  multivalued  func¬ 
tion,  and 

N(c-1(c(z) ) )  =  N(z ) , 

A  A 

where  z=z  or  any  depth  with  c(z)  =  c(z).  In  particular, 


and  equation  (4-1)  can  be  written  in  the  form  of  an  Abel 
integral  equation  (per  Reference  4-1): 


Notice  that  a  (z)  is  the  independent  variable 
(not  z),  and  that  the  equation  is  trivial  if  a  (z)  is 
constant.  In  fact,  assuming  continuity  of  c(z),  if  a  (z) 
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N(6; 


were  independent  of  z,  it  would  have  value  a(z)  =  0  and 
cos  a ( z )  =  1.  In  that  case,  equation  (4-1)  reduces  to 

tt/2 

i-i  r 

47  No  =  47  N^2  ^  ~  J  Ns<0;  2o^  cos0d0»  for  all  z, 

0 

and  any  directionality  function  Ns(6;z0),  suitably  normal¬ 
ized,  solves  the  integral  equation.  Again,  constant  sound 
speed  leads  to  a  degenerate  relationship. 


The  integral  equation  is  non-trivial  as  long 
as  a  ( z)  is  non-constant.  If  ot(z)  is  continuous  on  some 
interval  in  depth  and  if  N  and  c  are  suitably  smooth,  then 
the  Abel  equation  has  the  unique  solution  (in  the  class 
of  continuously  differentiable  functions;  Ref  4-1): 


zQ)  sin0  cos0 


For  v  =  cos  u, 

N  (6;Z  )  = 

s  '  o 


The  directionality  N  (9 ;z0)  can  thus  be  determined 
explicitly  from  the  depth-dependence  function  N(z)  via  the 
latter  formula.  Numerical  quadrature  and  differentiation 
are  usually  required.  There  are,  however,  special  forms 
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for  N(z)  in  terms  of  c(z)  which  allow  for  closed-form 
solutions.  Suppose,  for  example,  that  c(z)  is  not  constant 
and  that  N(z)  can  be  written  as: 


n 


N  (z)  =  A  +  A, 


c(zo)  ‘ 

X  4  •  t  X  A 

c(zo)' 

c(z) 

T  •  *  «  A 

n 

c(z) 

Then 


N  f0"1  [  v^~j)=Ao+  Alv  +  ...  +  An  v° 


and  thus 

N  (e  -Z  )  -  - - -  -4 

s  ’  2ir2sin0  cos0 


fll  V 


2d< 


50- v2 


2T:2sin0  cos0 


d0 


+1 


10=1 


_JL 

"  27 r2 


m=l 


+1 


E 

m=l 


(  2m  )  !  ttA  . 


(m-1) !  J2-: 


|n+l 

12 


See  Ref.  4-2 


E 

m=l 


2*m 

r 

( 2m+l ) ! 

i) 

2 (m-1) 

cos. 

2m 

2  *m 

A>m-1 

cos2™  0 


2m+l  _ 
cos  0 


(2m-l) ! 


cos203"1  0 


.  (4-2) 
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In  particular,  for  n=0  (i.e.,  depth-dependence  function 
is  constant),  the  noise  field  must  be  isotropic: 


Ns<0;V  -  4?  * 


Case  2  Asymmetry  in  9 


Suppose  Ns(0;zo)  is  not  symmetric.  Define 

ns(0;zo) 

g(e)  =  ¥TreTz~T  for  °-e$ir/2 


and 


V6:Zo}  '  2 


Nste;z0)-ms(-e;z0) 


iNs(e;zo) 


l+g(6) 


,  for  0  <  0<ir/2 , 


Then,  recalling  (2-7) 


N(z)  *  2tt 


r  -  / 

z)  J- 


-a(z) 

t/2 


d0Ns(0;zo)-F(0,z) 


/•tt/2  /*tt/2 

■  2*  /  ,  Ns(9;z  )F(0,z)d0+2ir/  g(  0  )NS(0  ;zo) -F(0  ,  z)d0  , 
«/a(z)  u  ./a(z)  &  ° 
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and  so, 


0(2)  -  4*  /  ”/2  N  (6,2  )F(6,z)d6. 

-Mz)  s  0 

A  solution  for  Ns(e;z0)  can  be  found  as  before 
in  (4-2).  Finally,  g(0)  can  be  estimated  by  tbe  sum  of 
the  volume  attenuation  along  the  additional  ray-path 
length  and  the  bottom  loss/attenuation,  if  any,  for  the 
additional  bottom  bounce  (see  Figure  4-1). 

4.2  Error  Minimization  Algorithms 

As  an  alternate  to  the  explicit  inverse  approach, 
consider  an  approximate  method  whereby  Ns(0  ;zQ)  is  assumed 
to  have  a  special  form,  depending  on  several  parameters. 
The  resulting  N(z)  calculated  from  (4-1)  is  compared  to  the 
actual  depth  function,  and  the  parameters  adjusted  to 
minimize  the  error.  In  the  illustration  below,  Ns(0;zo) 
is  assumed  to  depend  on  3  parameters,  but  the  approach  is 
directly  general izable  to  more  complicated  forms.  An 
intuitively  satisfying  observation  is  made  here:  the 

problem  is  well  defined  as  long  as  the  number  of  distinct 
samples  (i.e.,  those  at  distinct  sound  speeds)  of  N(z) 
is  at  least  as  great  as  the  number  of  unknown  parameters 
defining  Ns(0;zo). 


m 

Let  (Zi}i=i  be  a  sequence  of  depths  all  contained 
in  the  same  sound  channel  as  z0  with  c(zi)  >  c(z0)  for  all  i. 


I 
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Figure  4-1.  g(0)  is  approximated  by  additional 

attenuation  along  dashed  ray  path. 


Let  N(z^)  denote  the  omnidirectional  noise  level  at  depth 
z^.  As  was  shown  earlier,  if  the  noise  directionality  at 
zQ  is 


1Ni-  lel<es 
N2-  6s£ie^6B 

N3,  8b<|9|<»/2  , 


(4-3) 


then 


where 


Thus,  given  (NCz^)}  and  assuming  Ns  takes  form  (4-3), 
we  wish  to  find  non-negative  numbers  (N^,  N2,  N3)  such 
that  (4-4)  is  true  for  all  i. 
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First  note  that  when  only  three  depths  are  given, 
with  c(zi)  =£  c(zj)  for  i  *=  j ,  then  ignoring  any  error 
terms,  the  problem  is  simply  to  solve  a  system  of  three 
linear  equations  in  three  unknowns.  It  can  be  shown  that 
the  matrix  of  coefficients, 


M= 


°1 

ei-“i 

a2 

62-a2 

*-*2 

a3 

Va3 

*-*3 

is  non-singular  under  the  conditions  that  the  c(z^)  are 
distinct  and  that  cs  f  cb-  Hence  there  exists  a  unique 
solution,  namely 

(nA  /n(  Zx  )\ 

N2  )  =M-1  I  N(Z2)  I 
Ns/  \n(Z3)/ 

Now  when  more  than  three  depths  are  given,  the 
system  of  equations  is  overdetermined.  In  practical  appli¬ 
cations,  equality  in  (4-4)  for  all  i  is  virtually  impossible 
because  of  the  special  form  of  (4-3)  and  the  measurement 
accuracy  of  N(zj_).  Therefore,  some  sort  of  error  minimi¬ 
zation  process  must  be  used.  Two  such  processes  were 
considered:  (1)  a  least-squares  estimate  and  (2)  a  linear 
program.  Of  the  two  the  least-squares  estimate  requires 
fewer  computer  resources.  The  formulations  of  these  two 
problems  are  given  in  the  following  subsections. 
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4.2.1 


Least-Squares  Estimate 

Problem:  Find  the  least-squares  estimators  0^,02, 

03,  of  Ng,  for  function 


N(z) 


4tt 


Nx  ct  +  N2  (6-ct)  +  N3  (1-8) 


given  data  points  (aj.,  8i,N(zi),i  =1,  2,...,  m) 


Solution:  The  sum  of  the  squared  deviations  between 
the  measured  and  the  theoretical  values  of 
N( Zj[ )  is 


S  - 


Z)  I  N(2.)-4Tr  I  Nia±  +  N2  (B.-cu)  +  Ng  (X-e±)j | 


Then  setting  the  partial  derivatives 


as  as  as 
3N1'  aN2  ana  aN3 


to  zero,  we  obtain  nj,  n2»  and  n3  (  the  least-squares 
estimators  of  Nj,  N2  and  N3)  as  solutions  of  the  equations 

£N(Zi)a.=47r  N1£<x[  *  N£  £0^(8^)  +  N3  .£ai(  1-B± ) 

2^N(zi)(Bi“0ii  )=4tt  N£ai(ei-ai)-t  N2  E(ei"ai)2  +  N3  <E ( Bi~ai  )( 1_6i ) 
EN(zi)(l-8i)=47T  Nx  Ecx.d-B^  +  N2  £(  VV  a-B±)'+  Ng  £(1-B.)2 
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4.2.2 


Problem:  Define  yi  to  be  the  absolute  difference 
between  the  measured  noise  intensity  at  the  depth  and  the 
theoretical  noise  intensity  as  given  by  (4-4).  For  each  i 
assign  a  relative  cost  (penalty),  ci,  to  be  applied  to  the 
y^.  If  the  objective  is  to  minimize  the  total  cost,  then 
the  problem  can  be  solved  by  the  linear  program: 


Minimize : 


m 

i=l 


Subject  to: 

N(z±)  -  4tt  Nia.  +  n2(  6.-a.)  +  N3(l-6.  )  iyi>i=l, - ,m 

Nj  >  0  ,  j-1,2,3 
yi  >  O  ,  i-1,2, • •• ,m  , 

The  linear  program  is  advantageous  in  that  a  variety  of 
meaningful  options  for  the  cost  function  can  be  used 
(e.g.,  minimize:  max  Ci»yij).  Additionally,  constraints 
that  reflect  a  priori  knowledge  can  easily  be  included 
(e.g.,  N2  _>  2N3) . 
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Section  5 

EXAMPLES  OF  EFFECTS  OF  NOTCH 
FILLING  ON  DEPTH  DEPENDENCE 


In  subsection  3-2,  it  was  shown  that  when  the 
vertical  directionality  at  the  sound-speed  minimum  (c(zG)) 
takes  special  form 

h 

Ns(0;zo)  =jN2 

(N3 


,  o  < | e | <es 
,  es  <t«|<eB 

■  eB 


then  the  depth  dependence  function  is  given  by 


N(z)  =  4tt 


N„ 


sin0  +  No(sin0n  -  sinOl 

S  ZD  b 


N3(l-sin6B) 


where 


a 


arccos 

arccos 


c(z) 

c<zo> 

c(z) 

c(zo> 


cos0 

s 

COS0g 
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This  section  provides  examples  of  depth-dependence  functions 
in  a  variety  of  environments  for  vertical-directionality 
functions  as  defined  above.  Furthermore  these  examples 
illustrate  the  relative  effect  of  "notch  filling"  in  these 
environments . 

Figures  5-1  thru  5-4  are  for  a  typical  deep-water 
sound-speed  profile  with  a  well-defined  sound  channel.  At 
the  channel  axis  the  surface  and  bottom  grazing  angles  are 
approximately  10.3°  and  18°  respectively. 

In  Figure  5-lc  the  depth  dependence  varies  by 
about  2.5  dB^  reaching  a  minimum  at  the  axis  of  the  sound 
channel  and  a  maximum  at  the  surface  and  the  conjugate 
depth.  This  variation,  which  occurs  despite  the  assumption 
that  the  levels  N^,N2  and  N3  remain  constant  in  depth, 
is  the  effect  of  the  angle  transformation  given  by  (2-2). 
As  the  sound  speed  increases  away  from  the  axis  the  angles 
cs  and  eg  defining  the  vertical  directionality  are 
shifted  towards  zero  as  shown  in  Figure  5-ld.  The  results 
of  this  transformation  are  two-fold.  First,  the  dominant 
RSR  paths  are  given  more  weight  by  the  cos  0 -term  in  the 
calculation  of  omnidirectional  noise: 

N(z)  =  J*Ng(0;z)  cos0  d8 . 

Second,  the  angular  aperture  of  RSR  paths,  eB_es> 
increases . 


Figure  5-1.  Environment  One 

(a)  Vertical  Directionality  at  Channel  Axis 

(b)  Sound  Speed  Profile 

(c)  Depth  Dependence  Function 


Figure  5-1.  Environment  One 


(d)  Vertical  Directionality  Near  Critical  Depth 


Figure  5-2.  Environment  One 

(a)  Vertical  Directionality  at  Channel  Axis 

(b)  Sound  Speed  Profile 

(c)  Depth  Dependence  Function 


Figure  5-3.  Environment  One 

(a)  Vertical  Directionality  at  Channel  Axis 

(b)  Sound  Speed  Profile 

(c)  Depth  Dependence  Function 


Figure  5-4.  Environment  One 


(a)  Vertical  Directionality  at  Channel  Axis 

(b)  Sound  Speed  Profile. 

(c)  Depth  Dependence  Function 


For  this  particular  case,  at  the  channel  axis  the  RSR  aperture 
is  approximately  7.7°  while  at  the  critical  depth  it  is 
approximately  14.8°.  The  second  effect  is  generally  an 
order  of  magnitude  greater  than  the  first. 

Rather  dramatic  changes  in  the  depth  dependence 
occur  as  the  noise  level  in  the  notch  is  increased.  In 
Figure  5-2  where  the  notch  level  is  3  dB  down  from  the 
RSR  level,  the  depth-dependence  function  has  become  nearly 
constant.  As  the  notch  is  filled  further  (Figure  5-3),  the 
depth  dependence  becomes  concave  in  the  opposite  direction 
so  that  it  reaches  a  maximum  at  the  sound-channel  axis 
rather  than  a  minimum  as  in  Figure  5-1.  This  also  is  a 
result  of  Snell's  Law.  At  depths  where  the  sound  speed  is 
greater  than  at  the  axis,  some  of  the  shallower  arrival 
angles  in  the  notch  are  no  longer  seen.  In  Figure  5-1  the 
intensity  of  the  arrival  lost  is  negligible.  However,  as 
the  notch  is  filled  such  energy  becomes  significant. 

The  final  figure  for  this  environment,  Figure 
5-4,  demonstrates  the  damping  effect  on  the  depth-dependence 
function  caused  by  the  addition  of  more  energy  to  the 
bottom-bounce  paths.  This  is  as  expected  since  the  energy 
in  the  RSR  paths  is  then  a  smaller  percentage  of  the  total. 

The  second  environment,  Figures  5-5  thru  5-7, 
is  an  example  of  an  extreme  summer  profile,  with  a  strong 
thermocline  and  nearly  bottom-limited  conditions.  At  the 
channel  axis,  the  surface  and  bottom  grazing  angles  are 
approximately  15.4°  and  16.6°  respectively. 
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Figure  5-5.  Environment  Two 

(a)  Vertical  Directionality  at  Channel  Axis 

(b)  Sound  Speed  Profile 

(c)  Depth  Dependence  Function 


In  comparison  with  Figure  5-1,  the  depth  depend¬ 
ence  function  in  Figure  5-5  has  a  range  of  nearly  5  dB.  The 
increased  variation  is  caused  by  the  increase  in  the  notch 
width  and  the  narrowness  of  the  RSR  aperture  at  the  channel 
axis.  These  two  conditions  result  in  the  RSR  aperture 
varying  from  1.2°  at  the  channel  axis  to  6.3°  at 
critical  depth. 

Figures  5-6  and  5-7  again  demonstrate  that  as 
the  notch  is  filled  the  depth-dependence  function  completely 
changes  shape.  This  change  is  more  dramatic  than  in  the 
first  environment  because  proportionally  more  energy  is 
being  added. 

The  final  environment,  Figures  5-8  thru  5-10, 
further  demonstrates  the  importance  of  the  size  of  the 
RSR  aperture.  At  the  channel  axis  the  surface  grazing 
angle  is  7.3°,  smaller  than  either  of  the  first  two 
cases.  Yet  the  depth  dependence  function  in  Figure  5-8  has 
a  range  of  3  dB,  more  than  for  the  first  environment. 
This  happens  because  the  width  of  the  RSR  aperture  ranges 
from  1°  to  4°,  a  variation  which  is  proportionally 
greater  than  in  the  first  case.  The  effect  of  filling  the 
notch  in  this  case,  Figure  5-9  and  5-10,  is  similar  to  that 
of  the  previous  example. 
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Figure  5-8.  Environment  Three 

(a)  Vertical  Directionality  at  Channel  'Axis 

(b)  Sound  Speed  Profile 

(c)  Depth  Dependence  Function 


Figure  5-9.  Environment  Three 


(a)  Vertical  Directionality  at  Channel  Axis 

(b)  Sound  Speed  Profile 

(c)  Depth  Dependence  Function 


Figure  5-10.  Environment  Three 

(a)  Vertical  Directionality  at  Channel  Axis 

(b)  Sound  Speed  Profile 

(c)  Depth  Dependence  Function 


In  conclusion,  these  three  examples  suggest  that 
the  variability  in  the  omnidirectional  noise  as  a  function 
of  depth  is  dependent  on: 

(1)  the  width  of  the  horizontal  notch  in  the 
vertical  directionality  and  the  level  to 
which  it  is  filled, 

(2)  the  width  of  the  RSR  aperture  (particularly 
in  the  nearly  bottom-limited  cases), 

(3)  and,  to  a  lesser  extent,  the  amount  of  energy 
in  the  bottom-bounce  paths. 
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Section  6 

SUMMARY  AND  CONCLUSION 


Beginning  with  a  previously  documented  analyti¬ 
cal  result,  a  theoretical  relationship  between  ambient 
noise  vertical  directionality  at  a  single  depth  and  the  omni 
level  as  a  function  of  depth  has  been  derived  and  investi¬ 
gated.  This  earlier  result  was  an  approximation  for  noise 
vertical  directionality  as  a  function  of  depth  based  on  a 
simple  angle  transformation.  It  allows  omni  level  as  a 
function  of  depth  (restricted  to  depths  within  a  single 
"sound  channel  ")  to  be  estimated  in  terms  of  the  vertical 
directionality  at  a  single  depth  (the  channel  axis).  The 
principal  relationship,  and  those  that  follow  from  it,  are 
limited  by  assumptions  of  range-averaged  transmission  loss 
and  geometric  acoustics.  Extension  to  account  for  image 
interference  and  diffraction  have  been  discussed. 

Sections  3  and  5  described  the  specific  effects 
that  various  cononical  vertical  directionalities  have 
on  the  depth-dependence  function,  with  closed-form  expresr 
sions  derived  for  special  cases.  Section  5  further  investi¬ 
gates  the  effects  of  the  environment  for  one  of  the  canoni¬ 
cal  cases.  Of  particular  importance  is  the  relationship 
between  the  depth-dependence  function  and  the  amount  of 
energy  in  the'  notch  (the  null  near  the  horizontal  predicted 
by  range-independent  noise  models)  in  the  vertical  direc¬ 
tionality.  Examples  show  that  the  depth-dependence  can  be 
dramatically  affected  by  energy  added  to  the  notch  and  to  a 
slightly  lesser  extent  by  the  width  of  the  aperture  of  RSR 
arrivals. 
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The  inverse  relationship  (directionality  as  a 
function  of  omni  depth  dependence)  is  established  in 
Section  4  for  symmetric  directionality  functions.  A 
closed-form  expression  for  vertical  directionality  is 
found  for  depth-dependence  of  a  special  form,  and  two 
error-minimization  algorithms  for  approximating  the  vertical 
directionality  are  also  described.  Extension  to  the  asym¬ 
metric  case  requires  a  priori  knowledge  of  the  relationship 
in  vertical  directionality  between  positive  and  negative 
angles. 

As  discussed  in  Section  1,  the  motivation  for  this 
report  was  the  desire  to  use  existing  depth-dependence 
data  to  provide  a  clue  as  to  whether  a  notch  at  angles  near 
the  horizontal  occurred  in  the  vertical  directionality. 
Unfortunately  the  data  were  too  sparse  (seldom  more  than  two 
or  three  measurements  in  the  main  sound  channel)  for  the 
results  of  Section  4  to  yield  more  than  coarse  estimates. 
Search  for  better  data  and  a  definitive  evaluation  of  the 
results  are  indicated. 

Besides  evaluation,  future  work  should  address: 

•  a  priori  estimates  of  asymmetry  based  on 
bottom  interaction, 

•  error  minimization  algorithms  for  canonical 
directionalities  with  singularities  at  the 
surface  grazing  angles, 

•  the  impact  of  SII  and  wave  effects  on  noise 
depth  dependence, 

•  applications  to  transmission-loss  models. 
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